Limits on Preserving Quantum Coherence using Multi-Pulse Control 
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We explore the physical limits of pulsed dynamical decoupling methods for decoherence control 
as determined by finite timing resources. By focusing on a decohering qubit controlled by arbitrary 
sequences of 7r-pulses, we establish a non-perturbative quantitative upper bound to the achievable 
coherence for specified maximum pulsing rate and noise spectral bandwidth. We introduce numeri- 
cally optimized control "bandwidth-adapted" sequences that saturate the performance bound, and 
show how they outperform existing sequences in a realistic excitonic- qubit system where timing con- 
straints are significant. As a byproduct, our analysis reinforces the impossibility of fault-tolerance 
accuracy thresholds for generic open quantum systems under purely reversible error control. 
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Building on the discovery of spin-echo and multiple- 
pulse techniques in nuclear magnetic resonance [1], dy- 
namical decoupling (DD) methods for open quantum sys- 
tems [2] have become a versatile tool for decoherence con- 
trol in quantum engineering and fault-tolerant quantum 
computation. DD involves "open loop" (feedback-free) 
quantum control based on the application of a time- 
dependent Hamiltonian which, in the simplest setting, 
effects a pre-determined sequence of unitary operations 
(pulses) drawn from a basic repertoire. Physically, DD 
relies on the ability to access control time scales that are 
short relative to the correlation time scale of the inter- 
action to be removed. The reduction in decoherence is 
achieved perturbatively, by ensuring that sufficiently high 
orders of the error-inducing Hamiltonian are removed. 
Recently, a number of increasingly powerful pulsed DD 
schemes have been proposed and validated in the lab- 
oratory. Uhrig DD (UDD) sequences [3], for instance, 
perturbatively cancel pure dephasing in a single qubit 
up to an arbitrarily high order n while using a minimal 
number (n) of pulses, paving the way to further opti- 
mization for given sequence duration [4, 5] and/or spe- 
cific noise environments [6], to nearly-optimal protocols 
for generic single-qubit decoherence [7]. Experimentally, 
UDD has been employed to prolong coherence time in 
systems ranging from trapped ions [4, 5, 8] and atomic 
ensembles [9] to spin-based devices [10], and to enhance 
contrast in magnetic resonance imaging of tissue [11]. 

In a realistic DD setting, the achievable performance is 
inevitably influenced by errors due to limited control as 
well as deviations from the intended decoherence model. 
Since it is conceivable that both model uncertainty and 
pulse non-idealities can be largely removed by more ac- 
curate system identification and control design, some of 
these limitations may be regarded as non-fundamental 
in nature. Composite-pulse [12] and pulse-shaping [13] 
tcchinques can be used, for instance, to cancel to high 
accuracy the effects of both systematic control errors and 
finite-width corrections. We argue, however, that even in 
a situation where pulses may be assumed perfect and in- 
stantaneous, an ultimate constraint is implied by the fact 



that the rate at which control operations are effected is 
necessarily finite - as determined by a "minimum switch- 
ing time" T m i n for the available control modulation. Our 
goal in what follows is to rigorously quantify the perfor- 
mance limits to preserving coherence using DD as arising 
from the sole constraint of finite timing resources. 

We focus on the paradigmatic case of a single qubit un- 
dergoing pure dephasing due to either a quantum bosonic 
bath at equilibrium or classical (Gaussian) noise, and 
controlled through a sequence of instantaneous ir pulses. 
While representing an adequate idealization of realistic 
decoherence control settings [4, 5, 8, 9, 11], this problem 
is exactly solvable analytically [2, 3], enabling rigorous 
conclusions to be established. Our first result is a non- 
perturbative lower bound for the minimum decoherence 
error achievable by any DD sequence subject to a timing 
constraint r m i n , for noise spectra characterized by a finite 
spectral bandwidth uj c . Secondly, we show how to gener- 
ate "bandwidth-adapted" DD sequences that achieve op- 
timum performance over a desired storage time while re- 
specting the pulse-rate constraint, and demonstrate their 
advantages in a realistic excitonic qubit. Conceptually, 
our analysis highlights connections between DD theory 
and complex analysis of polynomials, and provides fur- 
ther insight into the fundamental capabilities and limi- 
tations of open-loop non-dissipative quantum control. 

Control setting. — Our target system is a single qubit 
whose dephasing dynamics in the quantum regime is de- 
scribed by a diagonal spin-boson Hamiltonian of the form 
H = H s <g> I B + H SB + Is® H B , with H s = uS z and 

H SB = S Z ® ^2(g k b k + g* k bl), H B = ^ u) k b\b k . 

k k 

Here, Is(B) denote the identity operator on the system 
(bath), S z = hcr z /2 is the spin operator along the quanti- 
zation axis, and b k (b k ) are canonical ladder operators for 
the fcth bosonic mode, characterized by a frequency cu k 
and coupling strength g k . If the bath is initially at ther- 
mal equilibrium at temperature l/(k B P), its influence on 
the qubit dynamics is encapsulated by the spectral den- 
sity function I(u>) = J^k Iff^l 2 ^^ — w fe)- Without loss 
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of generality, we shall assume that I{oj) decays to zero 
beyond a finite ultraviolet cutoff uj c . 

DD over an evolution interval [0, T] is achieved by 
applying a train of n instantaneous tt pulses (each im- 
plementing a Pauli o~ x operator) at times {tj}, where 
< ti < . . . < t n < T, and we also let to = and 
t n+ \ = T . While keeping the number of pulses n to a 
minimum may be desirable for various practical reasons, 
neither n nor the resulting sequence duration need to be 
constrained a priori. An arbitrary long duration T may, 
in fact, be needed for quantum memory. In contrast, infi- 
nite pulse repetition rates are both fundamentally impos- 
sible and undesirable as long as T > 0. Let the minimum 
switching time T m i n > lower-bound the smallest control 
time scale achievable by any sequence: 

r = mm(t j+ i - tj) > r min , j G {0, . . . , n). (1) 

If the system is initially in a nontrivial coherent su- 
perposition of S z eigcnstatcs, its purity in the presence 
of DD decays with a factor of exp(— 2x{t})> where the 
decoupling error X{t } > can be exactly expressed in 
the following form (see e.g. Eqs. (8c) and (10) in [3]): 

r°° f. 
X{t j} = / AH|/ ft} H| 2 d W , tj == — , (2) 
Jo T 

n 

f {ij} (u) = £(-l)V >T - e*» lUT ), (3) 
j=o 

and the "spectral measure" = 
2 coth(/3o;/2)/(a;)/cti 2 . In terms of the rescaled pulse 
times tj, Eq. (1) becomes tj+i — tj > 1. Physically, 
Eqs. (2)-(3) can also describe the purity decay resulting 
from pure dcphasing in the semi-classical limit, as due 
to stochastic fluctuations of the qubit energy splitting 
and experimentally investigated in [4, 5, 9]. In this 
case, H SB = and H s = [oj + £(t)]S z , where £(t) is a 
Gaussian random variable with a power spectrum S(ui) 
[14]. In order to evaluate X{tj}> it suffices to redefine 
X(u) = S(lu)/2ttlo 2 . The objective of DD is to minimize 
X{tj}- Our main problem then directly tics to the 
following: Given the fundamental constraint of Eq. (1), 
what is a lower bound on X{tj}? 

Non-perturbative performance bound. — A lower bound 
on X{tj] can be obtained by restricting the integral in 
Eq. (2) to a finite range [0, cj c ], with a tight bound ensu- 
ing if u c coincides with the spectral cutoff in either I(tv) 
or S(uj). We separate the dependencies of X{tj} upon 
the timings and the spectral measure A(u>) by applying 
Cauchy's inequality to the functions A 1 / 2 )/! and A -1 / 2 : 

Thus, the integral J^" y(u))\du), which is the Li-norm 
of the "filter function" fn.x over [0, u) c ], determines a 
worst-case lower bound on X{t } f° r °^ spectral densities 
A(w) for which the integral defining M{X} is finite. 



Interestingly, upon letting e tUT = z G C in Eq. (3) , the 
function takes the form of a complex "polyno- 

mial" P{t }( z ) with non- integer exponents. Such Miintz 
polynomials have been studied in the mathematical lit- 
erature, and a plethora of results (and conjectures) exist 
on their associated norm inequalities, zeroes, and multi- 
plicities [15]. The (now resolved) Littlewood conjecture 
in harmonic analysis [16] may be invoked, in particular, 
to lower-bound the Li-norm of f{t y 

X{t j} > ^C(logn) 2 , if u c t > 2tt, (5) 

with C = O(l). Also note that, regardless of oj c t, an 
upper bound follows immediately from Eq. (2): X{tj) < 
m {x} n2 , where jti{a} = J °° \(u))du>. Eq. (5) implies that 
in the "slow-control" regime where u> c t > 2ir, the DD er- 
ror worsens when more pulses are applied, and coherence 
may be best preserved by doing nothing. This reinforces 
how sufficiently fast modulation time scales are essential 
for achieving decohcrence reduction, as we discuss next. 

The "fast-control" regime (oj c t < 2ir) is implicit in 
perturbative DD treatments, where the filter function 
/{t }(w) is chosen to have a Taylor series that starts at 
(ut)™, so that X{t } remains small for sufficiently small 
values of ui c t. While this perturbative approach has been 
used for designing efficient DD schemes, it cannot lead 
to a lower bound on the attainable DD error in the pres- 
ence of a timing constraint. Consider for example UDD„ 
sequences, in which case tj = T sin 2 [irj/ (2n + 2)] for 
j = 1, • • • , n, and r = t\. If t is kept fixed, increasing 
n is only possible at the expense of lengthening the total 
duration as T(n) = 0[rn 2 ]. Irrespective of the fact that 
perturbatively the error scales as 0[(w c t)"], it carries a 
prefactor that grows too fast with n, eventually causing 
the perturbative description to break down [17, 18]. 

A non-perturbative lower bound may be established by 
directly mapping the Li-norm integral of f{ij] to the size 
of the corresponding Miintz polynomial P^ .y(z) over an 
arc of the unit circle of length uj c t. Theorem 2.2 in [19], 
in conjunction with Eq. (4), then implies: 

X{t ]} > ce~ a l^\ if uj c t < 2tt, (6) 

for some numeric constants c and a independent of r, w c , 
and {tj}. The bound in Eq. (6) is strictly positive for 
spectral measures of compact support. That it cannot be 
obtained by perturbative methods is manifest from the 
fact that it describes an essential singularity in lo c t. 

It is worth to further interpret Eq. (6) in the light of 
existing results. If the control rate is identified as the key 
resource that DD leverages for removing errors, a zero 
lower bound on X{tj} would allow, in principle, arbitrar- 
ily high DD accuracy to be achieved by using sufficiently 
long sequences with a fixed r m i n - that is, in analogy with 
fault-tolerant quantum computation [20], with a constant 
resource overhead relative to the noise-free case. Histor- 
ically, the impossibility of reliable computation with a 
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FIG. 1: (Color online) Decoupling error for UDD n sequences 
vs. cj c t, for a "flat" spectral measure A'°'(cj) = Q(u — w c ). 
The comparison curve denotes the general lower bound, Eq. 
(6), evaluated for a = 3, c = 1/2, chosen to approximate a fit. 



constant blow-up in resources (circuit depth) was estab- 
lished in [21] in the broader context of noisy reversible 
circuits, both classical and quantum. Therefore, our re- 
sults may be taken to reinforce the fundamental limita- 
tions of purely unitary quantum error correction, while 
explicitly characterizing the way in which such limiting 
performance depends upon the physical parameters. 

Achieving the performance bound. — Note that the 
minimum switching time r m j n enters Eq. (6) naturally, 
whereas both the total duration T and pulse number n 
are markedly absent from it. Thus interestingly, if the 
bound can be achieved, it should be possible to do so 
irrespective of how long T, provided that n is uncon- 
strained. We can show that the error associated with 
UDD sequences, Xn DD > saturates the fundamental limit 
in Eq. (6) in functional form although not necessarily 
in absolute sense (see also Fig. 1). This follows from 
noting that an upper bound to Xn DD m the presence of 
a hard spectral cutoff may be obtained from an upper 
bound to |/^ JDD (a;)|, by tailoring n to the bandwidth, 
n = uq rs l/(e 2 w c r) (see Remark 2.6 in [19]). This yields: 

< m {xy max |/„ U DI »| 2 < ^c'e^'^K 

where c' = 2/(ire 2 ), a' = 2/e 2 , and a similar functional 
form as in Eq. (6) is manifest. With t = t\ = r m j n fixed, 
the duration T of the "tailored UDD„" sequences scales 
as 0[l/(w 2 T min )], and the longest allowed r-value that re- 
sults in coherence improvement scales as l/(nu) c ). Thus, 
UDD provides no guarantee that the error reaches its 
absolute minimum and accessing the required r becomes 
increasingly harder as T grows. This motivates searching 
for DD sequences that can operate beyond the perturba- 
tive regime and retain their efficacy over the broadest 
range possible, up to l/{nuj c ) < r < l/w c - 

Various optimized DD strategies have been investi- 
gated for the qubit-dephasing setting under considera- 
tion. In "locally optimized" (LO) DD [4], optimal pulse 
timings are determined via direct minimization of the 
error X{tj} f° r a fixed target storage time T, whereas in 
"optimized noise filtration" (OF) DD, only the integral of 
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FIG. 2: (Color online) Decoupling error for BADD (dashed), 
LODD (dotted), and UDD (solid) sequences vs. total dura- 
tion T with the minimum interval r's indicated, for a dephas- 
ing exciton qubit operating at temperature 77K (see text). 
The search space for BADD and LODD covers up to n = 100 
pulses, whereas for UDD n < 20. See the Appendix for fur- 
ther detail and a comparison between different pulse timing 
patterns. 



the filter function is minimized [5] (see also [6] for a noise- 
adapted perturbative approach). While LODD/OFDD 
can access regimes where perturbative approaches are not 
efficient, they focus on matching the total sequence du- 
ration T as the fundamental constraint. However, this 
may fail to produce a satisfactory control solution if the 
timing constraint imposed by Eq. (1) is significant. 

To guarantee that such a fundamental limitation is 
obeyed, we introduce optimized bandwidth- adapted DD 
(BADD) sequences where both the minimum switching 
time and the total time are constrained from the out- 
set, see the Appendix for additional technical detail. 
We demonstrate the usefulness of BADD by focusing 
on the exciton qubit analyzed in [18], for which a spin- 
boson dephasing model with a supra-Ohmic spectral den- 
sity and a Gaussian cutoff is appropriate, I" c ' s (uj) = 
aw s exp(-w 2 /w 2 ) with s = 3, a « 1.14 x lCT 26 s 2 , 
uj c « 3radps _1 , and the need to avoid unwanted ex- 
citation of higher-energy levels enforces a timing con- 
straint Tmin « 0.1 ps [22]. The results are summarized in 
Fig. 2. Besides indicating the inadequacy of perturba- 
tive UDD for T > 1 ps, two main features emerge. First, 
as predicted by Eq. (6), the minimum error achievable 
by BADD is mainly dictated by r, largely independently 
of the total time T. Second, LODD performance is fairly 
sensitive to the timing constraint: for a fixed T (10 ps 
in Fig. 2), "softening" the constraint selects LODD se- 
quences that outperform BADD, the opposite behavior 
being seen if the constraint on the intended r is "hard- 
ened". Thus, a BADD protocol effectively optimizes over 
a set of LODD sequences where the timing constraint is 
only approximately met, consistent with intuition. 

In practice, an important question is whether the per- 
formance of a DD scheme is robust against uncertain- 
ties in the underlying spectral measure: in particular, 
sequences adapted to a presumed w c r m i n need not be ad- 
equate for the actual tv' c T m i n . Some illustrative results are 
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FIG. 3: (Color online) Purity loss, 1 - e 2x f'j>, vs. actual 
over presumed cutoff, oj' c /ui c , for the supra-Ohmic Gaussian 
spectral density (s = 3) corresponding to the exciton qubit. 
All sequences are adapted to T ~ lOps, r w O.lps. Varying 
the "actual" power law of the noise to s = 4 and s = 2 resulted 
in a qualitatively similar behavior (data not shown). 



depicted in Fig. 3 for sequences subject to the same tim- 
ing constraint, but applied to a setting where to' c =/= uj c . 
Clearly, a smaller cutoff oj' c leads to smaller decoherence, 
but much more so for perturbative UDD sequences. Ex- 
pectedly, the knowledge of the spectral density explicitly 
assumed in generating BADD and LODD results in far 
better coherence compared to OFDD and UDD, espe- 
cially when this knowledge is precise (lj' c /u) c = 1) or over- 
estimates the cutoff. Comparatively, BADD sequences 
appear to be more robust than LODD sequences when 
the cutoff is underestimated. 

Discussion. — Our mathematical description has re- 
lied on the solvability of the dephasing spin-boson model 
in the limit of instantaneous control pulses, however we 
expect similar timing-induced lower bounds to exist un- 
der more general conditions. In principle, non-Gaussian 
classical dephasing such as random telegraph noise could 
be addressed based on the exact solution presented in 
[23], whereas non-bosonic dephasing models of the form 
Hsb + Hb = S z (E)B z + Bo, could be tackled by matching 
the leading-order contributions in B z with the bosonic 
case studied here. Note, however, that bounded tim- 
ing resources do not prevent the DD accuracy bound to 
be zero in special cases - such as "monochromatic" or 
"non-dynamical" baths (Hb = 0), for both of which the 
length of the arc appearing in Eq. (6) vanishes. Similarly, 
"nilpotent" environments, where powers of the bath op- 
erators in Hsb and Hb vanish at some order, allow per- 
turbative DD schemes to achieve perfect decoupling, as 
perturbation theory becomes exact. For more "adversar- 
ial" environments, where Hsb is not restricted to but in- 
cludes single-axis decoherence, similar lower bounds must 
exist by inclusion. Elucidating the algebraic features re- 
sponsible for a finite vs. zero performance bound remains 
an interesting open problem with implications for quan- 
tum error correction in general. As opposed to pulsed 
control scenarios, continuous- time modulation subject to 
finite energy /bandwidth constraints has also been ex- 
plored for decoherence control [24]. Although, even for a 
purely dephasing qubit, finding the optimal modulation 



requires solving a non-linear integro-diffcrcntial equation, 
it would be interesting to quantify the extent to which 
the extra freedom afforded by continuous controls may 
improve the achievable performance lower bounds. 

We thank Michael Biercuk, Irene D'Amico, Daniel Li- 
dar, and John Prcskill for valuable input. Work sup- 
ported from the NSF through Grant No. PHY-0903727. 

Appendix: BADD Optimization Procedure 

The BADD optimization procedure uses the following 
input: the spectral measure A(w) used to characterize 
the environment (a continuous positive-valued function 
defined on positive real numbers); the minimum switch- 
ing time t used to constrain the pulse timings (a positive 
real number); the total duration T of the sequence (a 
positive real number larger than t). The output from 
BADD is the number of pulses n (a positive integer) and 
the pulse intervals {tj}™^ 1 (positive real numbers) that 
yield the smallest decoupling error while satisfying the 
following two constraints: 

n+l 

(1)t 3 >t, (2)^>=T. 

Both these constraints are linear in the variable 

x (") : = (ti,t 2 ,-- - ,t„ + i). 
Also notice that n is always bounded: 

1 < n < n max =T/t -I 

The quantity to be minimized is the coherence loss over 
time T, which is given by [see Eqs. (2)-(3) in the main 
text] : 

/>oo 

X (x) = / X(Lo)\Uoj)\ 2 (kj, 
Jo 

where 

n 

f x (u) = ^(-lVV tj "-e^ +1 "), 
i 

i=i 

The optimization problem described by x( x ) is highly 
non-linear and involves evaluation of x(x), itself an infi- 
nite numerical integral. Computationally, we replace the 
upper limit of infinity on the integral with a parameter 
Woo chosen such that 

/*QC /"Woo 

/ A(w)|/ X (w)| 2 dw« / X(uj)\f x (uj)\ 2 doj. 
Jo Jo 

In the soft-cutoff exciton qubit example considered in the 
main text, setting Woo = huj c is more than adequate. We 
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FIG. 4: (Color online) Pictorial comparison between the pulse 
timing patterns for the BADD and UDD sequences corre- 
sponding to T ~ 10, r = O.lps (points (a) and (b) in Fig. 
2 of the main text). Notice how the intervals of the BADD 
sequence (a) are compressed at the endpoints, but become 
effectively uniform mid-sequence. This resembles the (ana- 
lytical) interpolated DD protocol identified in T. E. Hodgson, 
L. Viola, and I. D'Amico, Phys. Rev. A 81, 062321 (2010). 



use the general-purpose routines of Matlab for evalu- 
ation and optimization of x( x )- I n particular, we hard 
code the form of A(w) into a Matlab function and eval- 
uate the integral for x( x ) using Matlab 's quadv func- 
tion (recursive adaptive Simpson quadrature). The nu- 
meric optimization is performed with the fmincov func- 
tion (linear constraints, nonlinear objective, and auto- 
matically chosen algorithm) with a numerical tolerance 
of 10~ 6 for both the optimal value and the constraint 
satisfaction. This procedure returns an optimal choice of 
intervals x^ n and a corresponding minimum error XmL 
for all possible values of n G {1, • • ■ , n max }. The mini- 
mum of {Xmili}"™" values identifies the BADD optimized 
solution. The numerical procedure for generating Fig. 2 
in the main text takes about 4 hours to finish using 3 
parallel jobs each running on a 2.7GHz quad-core cpu. 

In addition to the features discussed in the main text, 
the following qualitative features are worth emphasizing: 

• Theoretically, imposing constraints (l)-(2) together 
makes the search space compact even including the 
range of possible pulse numbers. Furthermore, Eq. 
(6) in the main text provides an explicit lower 
bound on the objective function, independent of 
the duration T as well as the pulse number. Both 
these facts provide a solid foundation for the BADD 



procedure. 

• The optimization always returns T\ = r, that is, 
the first interval uses the shortest available pulse 
interval. 

• When lu c t is sufficiently small, the overall BADD 
timing pattern approaches that of UDD pulse se- 
quences. Conversely, when ui c t is large, the tim- 
ing pattern approaches a periodic one, consisting 
of nearly equidistant pulses. This is depicted in 
Fig. 1 above for a BADD sequence (n = 40) and 
a UDD sequence (n = 15) operating at lo c t = 0.3, 
corresponding to the smallest allowed value for the 
excitonic qubit discussed in the main text. 

• Consequently, the optimal number of pulses n ap- 
proaches n max when tu c is large; for smaller tlu c 
the number of pulses required is smaller than n max . 

• The BADD optimized error does not vary signifi- 
cantly with the total time T, in line with the ex- 
pectation that a lower bound on decoupling error 
does not depend on T. 

The LODD optimization setting mentioned in the text 
differs from BADD since it uses n and T as the input 
(instead of r and T), and only constraint (2) is imposed 
in the minimization of the error x( x )- The OFDD opti- 
mization procedure is similar to LODD in this respect, 
however it uses a specific "flat" spectral measure in which 
X(lu) is equal to 1 for uj < u c and otherwise. That is, 
OFDD minimizes the integral of the filter function in- 
dependently of the actual spectral measure A, leading 
to a minimum error which always upper-bounds the one 
from LODD. As also discussed in the main text, for a 
given value of T the pulse sequences from BADD and 
LODD can effectively be matched by comparing them 
at (almost) equal r's. BADD sequences can thus be ap- 
proximated by searching LODD sequences with different 
pulse numbers for satisfying the minimum interval con- 
dition (1). 
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